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Abstract 

A novel non-perturbative gluon junction mechanism is introduced within 
the HIJING/B nuclear collision event generator to calculate baryon number 
transport and hyperon production in pA and AA collisions. This gluonic 
mechanism can account for the observed large mid-rapidity valence baryon 
yield in Pb + Pb at 160 AGeV and predicts high initial baryon densities at 
RHIG However, the highly enhanced A — A yield and the baryon transverse 
momentum flow observed in this reaction can only be partially described. 
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Recent data JT|-|5[ on p + A and A + B interactions at the CERN SPS has revealed 
a large degree of stopping and strange hyperon production in the heavy nuclear systems. 
The stopping is significantly under-predicted by models which assume that the primary 
mechanism for baryon transport is diquark-quark {qq — q) hadronic strings ||[7|. In this 
letter we implement a variant || of the baryon junction mechanism P,|T0| to address this 



problem. This implementation is now available as a Monte Carlo event generator, HIJING/B 
]TJ, a new version of the HIJING model [|T|. 

In Fig. 1, the rapidity distributions of valence protons and lambdas are compared to 
HIJING predictions. It is clear that even in p + S reactions, the mid-rapidity baryon distri- 
bution is largely under-predicted. In HIJING, as in the LUND Fritiof and dual parton 
model (DPM) ||14|| , the valence baryons cannot move inward from the beam and target frag- 
mentation regions by more than about 2 units of rapidity because of the assumed diquark 
fragmentation dynamics of the qq — q strings. 

Several models have considered additional mechanisms or string configurations to correct 



this problem. In the VENUS model [jig] , the stopping and strange baryon production is 
reproduced by adding an ad hoc double string mechanism whose parameters are adjusted 
to fit the data. The RQMD and UrQMD [[17]] event generators reproduce the observed 
stopping through incoherent multiple inelastic scattering of the valence (di) quarks. The 
difficulty with this idea is that the time dilation of inelastic processes at high energies leads 
to coherence which is absent in this model. 



In a new version of the DPM model [18], enhanced baryon stopping is more naturally 



achieved by introducing a diquark breaking mechanism flnj. In this approach, a gluon 
exchange mediates the break-up of the diquark by changing its color state from a {3} to a 



{6}. The cross section for this interaction is then estimated JTUJ to be 

^-(ph ^pX)^C {^jf \ ' cosh(y/2), (1) 

where C ~ 7 mb was estimated to fit the observed baryon stopping in pp collisions at the ISR 
[TIJ. The unusual cosh(y/2) rapidity dependence and energy dependence follow here 



from the Regge motivated 1/y 7 ^ distribution of the valence quarks in the proton at small 
x. The resulting configuration is a q — q — q string where the baryon number is associated 
with the middle quark. In nucleus-nucleus collisions, diquarks are assumed to remain broken 
once they transform into a {6} color state. Thus, in A + B collisions, the probability for 
diquark breaking should be P^b = 1 — (1 — PdbYi where n is the number of collisions 
and is the diquark breakup probability in a nucleon-nucleon collision. For a heavy 

nucleus, very few of the diquarks survive. The addition of this diquark breaking component 



was shown in [18| to reproduce the observed stopping in AA collisions at the SPS. However, 



this mechanism alone could not account for the observed strangeness enhancement without 



including additional final state interaction effects |19| . 

In this letter, we consider another Regge motivated mechanism || which also leads to 
eq. ([!]) but which differs from the diquark breaking model in multiplicity and strangeness 
enhancements. This mechanism is motivated from the non-perturbative gluon field con- 
figuration (the baryon junction) that appears when writing the simplest gauge invariant 
operator for the baryon in SU C (3); 
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Here, the baryon junction is the vertex at xj where the three gluon Wilson lines link the 
three valence quarks to form the gauge invariant non-local operator. In a highly excited 
baryonic state, the Wilson lines represent color flux tubes. When these strings fragment 
via qq production, the resulting baryon will be composed of the three sea quarks which are 
linked to the junction while the original valence quarks will emerge as constituents of three 
leading mesons. Being a gluonic configuration, it was proposed || that the junction could 
be more easily transported into the mid-rapidity region in hadronic interactions. 

In Regge phenomenology, allowing for the possibility of baryon junction exchanges in pp 
and pp scattering is taken into account by adding new Regge trajectories (M J ). Here, the 
exchange of the leading M$ trajectory is characterized by a three jet event in contrast to 
the two jet event resulting from the Pomeron exchange. Phenomenologically, the difference 
between the pp and pp topological cross sections, Act = cr(pp) — cr(pp), is given by the 
contributions of the p, uj and Mq Regge exchanges || . Analysis of the energy dependence 
and the multiplicity distributions of pp and pp data shows that Act oc s -1 / 2 [^C§ and that the 



difference, Act, may be largely associated with annihilations [pTJ] . Upon associating Act with 
the annihilations or the M J exchanges, the intercept of the leading term, a M j(0), can be 

estimated from the energy dependence of Act, i.e. Act oc s -1 / 2 ~ s M ° ^ , giving a M j(0) ~ 
1/2. This value of a M j(0) ~ 1/2 is also consistent with estimates of the intercept using 
multi-peripheral model approximations [2"2,E| . In addition, various models which associated 



Act with three sheet, Mq type annihilation events were able to explain the observed mean 
multiplicity, n BBann ~ \n BBscatt M, the observed multiplicity distributions of Act 



and the experimental values of the ratio R n = (cr(pp) n — (y{pp) n ) / a{pp) n , where n is the 
prong number ||20|1 . Recently, Kharzeev || used the exchange of the Mq Reggeon in baryon 
production to propose a new baryon stopping mechanism. This mechanism is able account 
for the excess of mid-rapidity valence baryons in pp at ISR energies and was suggested to 
provide substantial stopping at RHIC and the LHC In this letter, we implement this new 
baryon stopping mechanism into the Monte Carlo event generator HIJING/B, and tests its 
implications for RHIC. 

As was shown by Kharzeev || , the energy and rapidity dependence of the inclusive baryon 
production at mid-rapidity can be obtained using Mueller's generalized optical theorem WA 



in the double Regge limit. Here, the exchanges of a Pomeron and a Mq Reggion lead to the 
following form for single mid-rapidity baryon production, 

d 3 n-(V / <s \V 4 

Eb=£— = C B f B (m*) (^) cosh(y/2) . (3) 



d 6 p B \ S 

where C B is a constant that reflects the couplings of the Reggeon and Pomeron to the proton, 
/^(m 2 ) is an unknown function of m t and sq ~ 1 GeV is a Regge energy scale. The cosh(y/2) 
rapidity dependence and the energy dependence follow from the assumed intercept, 

a M j(0) ~ 1/2. In contrast to the diquark breaking model in the DPM, the multiplicity || 
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of these events is enhanced by a factor of 5/4 while the strangeness content is enhanced by 
a factor of 3. The enhancement factor of 3 of the strangeness content of the baryon allows 
for the unique possibility of producing S = — 3 Q~ baryons. 

We note that the value of the junction intercept, a M .i(0) = 1/2, has been criticized by 



the authors of the diquark breaking mechanism ||25|| . Following the phenomenology which 
motivated the diquark breaking mechanism, they show that the exchange of two gluons in 
the {10} color state (Decameron) leads to baryon number transport with baryon junctions. 
However, in contrast to the above model using the exchange of M J Reggeons, their model 
leads to a constant, energy independent annihilation cross section []2"5| , o"{i } — 1 — 2 mb, and 
a uniform, rapidity independent, inclusive cross section for mid-rapidity baryon production 
|10| , da(pp — > pX)/dy ~ 0.1 mb. They also fit |25|] the multiplicity distributions of Act at 
Ei a b = 10, 20, 50 and 100 GeV to argue that Aa is dominated by the p and u exchanges, 
while the three sheet annihilation component only contributes a small energy independent 
cross section of 1 — 2 mb. The energy dependence of this small three sheet annihilation cross 
section would then imply that a M j(0) ~ 1. 

In HIJING/B fll|] , we implement Kharzeev's || model through a "Y" string configuration 
for the excited baryon. For reactions without junction exchange, standard qq — q strings are 
used. The baryon is resolved around the junction via qq production and the resulting three 
beam jets are fragmented as q — q strings. The hard processes are modeled as kinks in one 
of the three qq strings. A value of obj = 18 mb is taken to reproduce valence proton data 
26| from p + p collisions at 400 GeV/c incident momentum (^ l /so = 27 A GeV). In order 



for the three beam jets to have sufficient phase space to decay, the junction exchange is 
only allowed if the invariant mass of the excited "Y" configuration exceeds m > 5 GeV. At 
SPS energies, this kinematic constraint considerably limits the number of junction exchanges 
allowed, reducing its effective cross section to ~ 9 mb. Like the modified DPM, we assume for 
multiple collisions that the baryon junction remains stopped in subsequent soft interactions 
once it has been exchanged. The of the junction baryon is obtained by adding the pr 
of the three sea quarks. We note that the present version of HIJING/B does not include 
final state interactions, as our main interest is to test the extent to which initial state 
non-equilibrium dynamics can account for the observed strangeness and p± enhancements. 

In Fig 1, the valence proton rapidity (dN/dy p — dN/dyp) and valence hyperon rapidity 
(dN/dyx — dN/dyx) distributions of HIJING and HIJING/B are compared with minimum 
bias p + S data at 200 AGeV and central Pb + Pb (published and preliminary) data 
at 158 AGeV While HIJING severely under-predicts both the observed stopping and 

the hyperon production, the junction physics of HIJING/B sufficiently enhances both the 
baryon stopping and the hyperon production at mid-rapidity, reproducing the p + S results. 
However, as seen in Fig Id, even with the sizeable strangeness enhancement factor of 3, this 
gluonic mechanism can only partially account for the observed hyperon enhancement in the 
heavy systems such as Pb + Pb. 

The pt distributions of the protons are calculated in HIJING/B for p + Pb at beam 
momenta of 450 GeV/c and S + Pb at beam momenta of 200 GeV/c in the rapidity interval, 
2.3 < y < 2.9. The itlt distributions were then fit for itit — m p < 0.68 GeV with the 
functional form Ed 3 N/dp 3 = Aexp(— ttlt/T), where A is a normalization constant and T is 
the inverse slope. The inverse slopes were found to be T w 151 MeV for p + Pb and T m 168 
MeV for S + Pb. The measured values from the NA44 collaboration [Hf are T = 195 ± 5 
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MeV for p + Pb and T = 256 ± 4 MeV for S + Pb. Although the < p\ > of the junction 
baryons is enhanced by the a factor of 3 over the normal qq — q string, the contribution 
from this mechanism alone is insufficient to reproduce the observed enhanced flow. This 
results from the competition between the different pt distributions which arise from the 
different production mechanisms. If one-half of the protons emerge from qq — q strings 
where Td q ~ 130 MeV, while the other one-half come from junction configurations where 
T BJ = V3T dq » 225 MeV, the effective inverse slope is only T e f / ~ 1 73 MeV in the measured 
region. The inability of this rather strong initial state non-equilibrium dynamical mechanism 
along with the already included Cronin effect to account for the observed transverse baryon 
flow provides evidence for its possible origin as due to final state interactions. We note also 
that the above junction dynamics alone provides no mechanism to account for the observed 
enhancement of the anti-proton px • 

The impact parameter dependence of this stopping mechanism is studied in Fig 2 for 
the net valence baryons, B — B = (p — p) + (n — n) + (A — A) . For the impact parameters 
6 = — 3 fm, b = 4 — 5 fm and 6 = 7 — 8 fm, a strong to moderate degree of baryon stopping 
is observed. However, at b = 10 — 11 fm, the degree of stopping has decreased and the shape 
suggests semi-transparency. Recent data |27j has shown the suppression of J/ty for impact 
parameter of b < 8 fm. Measurements of the impact parameter (E?) dependence of baryon 
stopping would be of interest to test if the anomalous J/ip suppression |27j threshold at 
b ~ 8 fm is correlated with the onset of greater baryon stopping. In the baryon junction 
exchange picture, a large degree of stopping is directly correlated with an enhanced gluonic 
field intensity at mid-rapidity that could partially be the cause of the ionization of c — c 
pairs. 

The predictions of this model for the valence proton and lambda rapidity distributions 
in Au + Au collisions at RHIC energies (y/s = 200 GeV) are shown in Fig 3. HIJING/B 
predicts approximately twice the initial number of valence protons and five times the initial 
number of valence hyperons of HIJING at mid-rapidity leading to a prediction of twice the 
initial baryon density, p(t ) ~ 2p ~ 0.3/fm 3 . Previous predictions for RHIC assuming 
idealized zero baryon chemical potential scenarios should therefore be re-examined. 
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FIG. 1. HIJING (solid) and HIJING/B (dashed) calculations of the valence proton and hyperon 
rapidity distributions are shown for minimum biasp+5 collisions at 200 AGeV and central Pb+Pb 
collisions at 160 AGeV. The data are from measurements made by the NA35 NA44 and 

NA49 I0I collaborations. 
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FIG. 2. Calculations of the net valence baryon rapidity distributions using HIJING/B are given 
for the impact parameter windows of b = — 3 fm, 5 = 4 — 5 fm, b = 7 — 8 fm and b = 10 — 11 fm. 
In this calculation, the net valence baryons are defined as B — B = (p — p) + (n — n) + (A — A). 
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FIG. 3. Predictions for the initial valence proton rapidity distribution (upper two curves) and 
for the initial valence hyperon rapidity distribution (lower two curves) are given for Au+Au colli- 
sions at E cm = 200 AGeV by HIJING (solid) and HIJING/B (dashed). 
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